d1<-read.table("SureSelect_All_Exon_V4+UTRs_with_annotation_hg19_chr22_q13_2.bed",sep="\t")
d2<-read.table("SureSelect_All_Exon_V4+UTRs_with_annotation_hg19_chr22_q13_2_modified.bed",sep="\t")
100000*sum(d2[3]-d2[2])/sum(d1[3]-d1[2])
savehistory("script_callength.R")
